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1. Introduction 

Simulations of lattice QCD in the so called e regime [JTJ| of the chiral expansion allow in princi- 
ple the extraction of physical parameters, like decay constants and electroweak effective couplings. 
Simulations in the e regime are not exclusive to lattice actions with exact lattice chiral symmetry. 
While topology certainly plays an important role in this extreme regime, here we aim at simulations 
which sample all the topological sectors. The motivation for this is that the original expressions are 
given for the situation where topology is summed over. 

To fix the notation we recall that we use as the lattice action for Nf = 2 degenerate flavours 

s[x,x,u}=s G [u}+s F [x,x,u], (i-i) 

where So denotes the tree-level Symanzik improved gauge action [0] 



S G [U] = ^-£< b £ MeTr [l-p( lxl )(x ;i u,v)j +b { £ ReTr [l -p( Zxl \x,H,v) 



(1.2) 



Dw = ^{7p(V M + V*) - fl V* VjJ + zno, (1.4) 



with the normalization condition bo = 1 — 8Z?i, and b\ = — X. The fermionic part of the action is 
given by the Wtm action [^, ||] 

Sf\X,X,U] =a 4 £z(x)[D w + W 5 t 3 ]*M, (1-3) 

X 

where 

" ^IfM^-r V ^~ UV M V MJ 

, V* are the standard gauge covariant forward and backward derivatives, mo and /i q are respec- 
tively the bare untwisted and twisted quark masses (for a recent review on Wtm and unexplained 
notations see ref. |^]). 

The chiral phase diagram of Wilson-like lattice actions, due to lattice artefacts, is different 
from the one in the continuum [Q, [7|, §]. In particular for our choice of the gauge and fermion 
action, the so called first order Sharpe-Singleton scenario [^] takes place. In this case at fixed 
lattice spacing a, if we set wipcac = where 

I x (^oAg(x)P a (0)) 
fflpCAC= 2I^(f a = XX (L5) 

the values of the twisted mass which can be simulated, are bounded from below by a value \i z 
which is proportional to fl ,2 Ag CD [10, 1 1 , 12]. This result is obtained if one analyzes the potential 



of the chiral lagrangian describing the low energy properties of the theory taking into account also 
the cutoff effects of the lattice action up to 0(a 2 ). In particular this analysis is performed in infinite 
volume. 

It is well known ^ that in continuum QCD if one sends the quark mass to zero keeping the 
size of the volume fixed, chiral symmetry is restored and no phase transition appears, i.e. the 
dependence of the chiral condensate on the quark mass is smooth. It is plausible to expect that 
at finite lattice spacing the same mechanism takes place, and in particular the critical point /i c 
smooths out. One way to show this is to include the effects of the non vanishing lattice spacing 
in the analysis of [jjj, and to study the mass dependence of the chiral condensate in the e regime 
including the 0(a 2 ) effects. This work is currently in progress. 
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2. Algorithm 

To simulate dynamical Wtm Nf = 2 degenerate quarks in the e regime, we propose to use a 



PHMC algorithm [ |13| , |14] , |15| , |16Q . We refer to these references for more details on the algorithm, 
and here we just shortly summarize the main ingredients. The expectation value of a generic 
observable & is given by 

= ^ / ®P\^ Sa[V] det (<2<2 f [U]) &[U] Q = 75 [Dw + i/itf] , (2.1) 

with Q being a single flavour operator, which in the following is normalized to have the biggest 
eigenvalue equal to one. It can be computed splitting up the determinant in two terms 

where P„ g (QQ?) is a polynomial approximation of order n for [<2<2^1 1 restricted in the region 
of the eigenvalue spectrum [ft, 1]. The main idea is to split the eigenvalue spectrum of the Wtm 
operator into parts which are then treated by either incorporating them in the update step of a 
simulation algorithm or by taking them into account in a reweighting procedure. Following [Q|, 
and leaving out the stochastic estimate of the correction factor, the generic expectation value can 
be written as 



, W = det [QQ*P n ,e (ee 1 ")] - EI PAe &)] . (2-3) 

\ W ' P Xi<e 

and the reweighting factor W is computed exactly evaluating N% ~ 20 smallest eigenvalues of 
the Wtm operator where A; denotes an eigenvalue of QQ} . This strategy should allow a better 
sampling of the configuration space, when the eigenvalues of the Wtm are particularly small. To be 
specific in the exploratory runs we have performed we have used a polynomial approximation with 
n = 380 and a cutoff in the spectrum e = 5 • 10~ 5 . The bare parameters used in the simulations are 
J8 = 3.9, ajUq = 5 • 10~ 4 and k = 0.160856 @. In fig. | we summarize first preliminary results 
concerning the quality of the runs. In the first stripe we plot the MC history and distribution of 
the reweighting factor W. We observe a very smooth behaviour and a well behaved distribution. 
This has to be compared with the MC history and distribution of the smallest eigenvalue of QQ} , 
which is very small, but with fluctuations which are under control. In the remaining three stripes 
we compare the MC histories and distributions of the pseudoscalar density correlator ( p~T| ) at time 
slice xo = T/4 for the PHMC without and with reweighting factor and for the mass preconditioned 
HMC with multiple time scale integrator jl8|] (mt-mHMC ) . This plot shows the different sampling 
of the configuration space by the two algorithms, and in particular the HMC seems to suppress 
contributions coming from the very low eigenvalues.These plots also show the expected advantages 
of the PHMC algorithm: better statistics for the large fluctuations and supression of exceptional 
fluctuations by the exact reweighting factor. 

3. Sampling 

To check that the algorithm correctly samples all the configuration space we have measured 
the topological charge v using the field theoretical definition (proportional to FF) on cooled and 
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Figure 1: In the first stripe we plot the MC history and distribution of the reweighting factor W. In the 
second stripe we plot the MC history and distribution of the smallest eigenvalue Xq of QQ' . In the remaining 
three stripes we compare the MC histories and distributions of the pseudoscalar density correlator at time 
slice xq = T /A for the PHMC without and with reweighting factor and for the mt-mHMC. 



APE or HYP [20] smeared gauge configurations. We emphasize that the aim of this exercise 
is just to have a first impression on how the algorithm is sampling the configuration space and it 
is not intended to be an attempt to compute physical quantities like the topological susceptibility. 
If we consider the sample of gauge configurations generated by the polynomial P„g we obtain 
distributions of the topological charge as the ones in fig |2[ In fig. |3| we compare the MC histories 
and distribution of the topological charge and of the minimal eigenvalue of QQ ' obtained with and 
without reweighting. We observe that the topological charge distribution N v , after introducing the 
reweighting factor, has a qualitatively different shape. 



4. Results 



A first preliminary result is given by the determination of the pseudoscalar density correlator 
Cp(x ) = -a^{P a (^)P a {0)) a = 1,2, (4.1) 
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Figure 2: Distribution of the topological charge obtained with the field theoretical defintion on cooled and 
HYP smeared gauge configurations. 
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Figure 3: Comparison of the MC histories and distributions of the topological charge, obtained with one 
level of HYP smearing. We observe that the topological charge distribution N v , after introducing the 
reweighting factor, has a qualitatively different shape, even if it shows that the algorithm is able to sam- 
ple many topological sectors without any problem. 



where 



P cl {x)=%{x)y 5 -x{x). (4.2) 
In the £ expansion we have a relation between the chiral condensate and Cp{xq) |^l|, ^2|]. The 



equivalent relation at fixed topology can be found in [|23|]. The correlation function is expected to 
be constant up to a curvature in the Euclidean time given by the higher order corrections. In fig. [| 
we show the Euclidean time dependence of the correlation Cp(xo), with the fit result and range. 
Inserting the values of the lattice spacing and renormalization factors 

a = 0.0855(5)(4), zJF(jit = 2GeV) = 0.46(1)(2) (4.3) 
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Figure 4: Dependence on xo/a of Cp(xo) in the e regime using Wtm fermions. A parabolic fit and the range 
used is also shown. 



computed by the ETMC [|24j, |2j, |26p we obtain for the chiral condensate the following determina- 
tion 

= 2GeV)\ = (264(12)(4)+2°MeV) 3 (4.4) 

where the first error is the summation in quadrature of all the statistical errors involved in the 
determination, the second error is an estimate of the systematic error coming from the error in the 
determination of Zp, and the third is an estimate of the higher order corrections coming from the e 
expansion. 



This value compares nicely with the p regime determination of ETMC [27] and of JLQCD [28] 
and other previous determinations summarized in tab.2 of ref. [E5]. 



5. Outlooks 

In this proceeding contribution we have shown that with Wtm fermions, simulations in the e 
regime are feasible. It would be interesting to see if in the current simulations we have a physical 
volume which is big enough to perform a safe matching with chiral perturbation theory. Work 
on improving the basic PHMC algorithm we have used so far is in progress. Improvements on 
the quality of the correlation functions could be achieved with stochastic sources and low mode 
averaging. Interesting further investigations concern, e.g. the distribution of the topological charge, 
and the interplay between the chirally breaking 0(a 2 ) effects and the contribution of the zero modes 
to physical observables. Investigations in these directions are in progress. 
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